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Abstract 

We consider the problem of stable matching with dynamic preference lists. At each time- 
step, the preference list of some player may change by swapping random adjacent members. 
The goal of a central agency (algorithm) is to maintain an approximately stable matching, 
in terms of number of blocking pairs, at all time-steps. The changes in the preference lists 
are not reported to the algorithm, but must instead be probed explicitly. We design an 
algorithm that in expectation and with high probability maintains a matching that has at 
most 0((log'n)^) blocking pairs. 


1 Introduction 

In the world of massive and distributed data, it is hardly reasonable to assume that data are 
static. Yet, one must design algorithms that maintain a solution for a given problem that is 
(approximately) consistent with the requirements, e.g., a permutation that is almost sorted. 
Thus, it is important to design algorithms and data structures that are robust to changes in 
their input, i.e., they produce an output with some performance guarantee (quickly). 

There are a few different dynamic data models that have been considered. The area of 
dynamic graph algorithms consists of maintaining some property or structure, such as connec¬ 
tivity, matchings, or spanning trees, even when the underlying graphs are changing |EGI991 
lURlOl INS131 [GP13j . Here, it is assumed that the changes to the graph may be arbitrary, but 
are reported to the algorithm; and the focus is on designing data structures and algorithms that 
adapt efficiently (typically in terms of computational time) to changes in the input. The area of 
streaming algorithms studies the setting where the data can only be accessed as a stream and 
the focus is on producing the desired output with highly space-efficient procedures (typically 
poly-logarithmic in the size of the input). In the area of online algorithms, one must design 
procedures that, even when data is revealed bit by bit, produce an output that is competitive 
with algorithms that see the entire input at once. 

Recently, Anagnostopoulos et al. [AKMUlT] proposed the evolving data model to take into 
account the dynamic aspects of massive data. In this model, the changes to the data are not 
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revealed to the algorithm; instead, an algorithm has query access to the data. However, it 
is assnmed that the changes to the data are stochastic, not adversarial. In this setting, the 
focns is not on compntational complexity (which is allowed to be polynomial at each time- 
step), bnt qnery complexity, the nnmber of probes made by the algorithm. Anagnostoponlos 
et al. |AKMUfT] stndied the problem of maintaining the minimum element of a permntation 
and an approximately sorted permntation, motivated by questions such as maintaining high 
(page)-ranked pages. In their setup, a permutation evolves by choosing a random element and 
swapping it with an adjacent element. In later work, Anagnostoponlos et al. |AKM^12 studied 
evolving graph models and problems such as s-t connectivity and minimum spanning tree. 

In this work, our aim is to bring this notion to game theory starting from the basic problem 
of computing a stable matching. In other words, we introduce the notion of evolving agents, 
who may not report any updates to their strategy (or preferences) without an explicit request. 
In the stable matching problem, the input consists of two sets A, B of equal size, and for each 
member a total order (preference) over members in the other set. Given a matching between A 
and B, a pair (x, y) with x € A and y a B is blocking if they prefer each other to their matches. 
A matching is stable if there are no blocking pairs. Gale and Shapley showed that a stable 
matching always exists and can be found by an efficient algorithm [GS62j . We consider the 
setting where the preference lists evolve over time. The preference lists can evolve over time, 
by swapping adjacent elements. More precisely, while the algorithm can perform one query per 
time-step, we assume that a total number of a swaps events also occur, where a = 0(1) is 
some hxed parameter, called the evolution rate. This assumption is similar to previous works 
and models the critical regime: with less evolution events the input is basically static, and with 
more the input evolves too fast. The goal is to maintain a matching that has few blocking pairs. 

We summarize our results as follows. All three statements hold in expectation and with 
high probability. 


(i) Using the results of Anagnostoponlos et al. |AKMU iT] for sorting permutations, we 
design an algorithm that maintains a matching with at most 0(nlogn) blocking pairs, 
at all time-steps after roughly the first n^logn steps iTheorem 13.11) . Also, we observe 
that any analysis that uses their method as a black-box, cannot improve on this bound 
(Remark 13.ip 0 

(ii) In a restricted setting, where only one side, say the B side, has evolving preference lists, 
and if the A side has uniform random permutations as preference lists (known to the central 
agency), we design an algorithm that maintains a matching with O(logn) blocking pairs 
at all time-steps after roughly the first nlogn steps (Theorem [32]) • 

(iii) Finally, we design an algorithm in the general setting, that maintains a matching with 
at most 0((logn)^) blocking pairs at all time-steps after roughly the first n^logn steps 
(Theorem 14.ip . 


2 Preliminaries 

In the rest of the paper, n > 1 denotes an integer parameter and [n] the set of integers 
{l,2,...,n}. For a non-negative random variable X, parametrized by some integer n, we 
write “A = 0(/(n)) in expectation and with high probability^^ when for any constant c there 
exist constants no,c',c" > 0 such that E[A] < c'f{n) and P(A > c"/(n)) < for every 
integer n > uq. 

^We don’t rule out the possibility that a more fine-grained analysis of the algorithm may give better bounds; 
instead we design new algorithms. 
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2.1 Stable Matching 

We only consider the bipartite stable matching problem, also known as stable marriage. There 
are two sets of players A and B, with |A| = \B\ = n. Each player x G A {y ^ B) holds a 
permutation of B (A), denoted tIx '■ B ^ [n] (vr^ : A —)■ [n]) indicating their preferences over 
players in the set B (A). Thus, for y G B, TTxiy) denotes the rank of y in x’s preference list 
(where 1 is the highest rank). 

Let M : A ^ B he a matching (a bijection). A pair (x, y) is said to be blocking if y M{x) 
and X -<T^y M~^{y), where 2 ; z' indicates that 2 ; is ranked higher than z' according to 
permutation vr [i.c., 7r(2;) < 7r(2;')). In words, x prefers y to M{x) and y prefers x to M~^{y). 

A matching M is stable if there are no blocking pairs. Then the stable matching problem 
is to find a stable matching given preference lists {vr^ : 2 : G A U B)}. Gale and Shapley |GS62j 
proved that a stable matching always exists, and gave an algorithm that given the preferences 
lists as input finds a stable matching in O(n^) time. 

The Gale-Shapley algorithm is simple to describe. Only players in the set A make proposals. 
Initially all players are unmatched. Let M denote a partial matching at some point. If there 
is an unmatched player x G A, x makes a proposal to y G B, where y is the highest-ranked 
player in tTx to whom x has not yet proposed. If y is unmatched, or prefers x to M~^(y), then 
y accepts the proposal and we set M(x) = y. In the latter case, the agent previously matched 
to y, i.e., M~^{y) before M was updated, becomes unmatched once more. Gale and Shapley 
showed that this algorithm always results in a stable matching. 

Wilson |Wil72j studied the problem where all the preference lists are independent and uni¬ 
formly random permutations; in this case, he showed that the number of proposals made by 
the Gale-Shapley algorithm is 0(n log re) in expectation and with high probability (see also 
|MR95j L In fact, only the proposing side needs to be random in their statement. We provide 
a proof sketch for completeness. 

Theorem 2.1 f [Wil72] l. If the permutations {tTx '■ x G A} are ehosen randomly, the number of 
proposals made in the Gale-Shapley algorithm (where only A makes proposals) is O(relogre) in 
expeetation and with high probability. 

Proof sketch. Following the proof in |MR95j (see also |Knu97j ), analyze an alternative procedure 
where every proposal is uniform over the whole of B. If it happens that x G A proposes to a 
y G B that has already rejected x, then a rejection is guaranteed. It is not hard to show that the 
number of proposals such an algorithm makes stochastically dominates the number of proposals 
of the classical algorithm. Next, by the method of deferred decisions, fix the randomness in 
the algorithm only when needed. Then observe that the number of proposals is equal to the 
number of coupons chosen in the coupon-collector’s problem. □ 

2.2 Model for evolving input 

A general framework for studying dynamic data was introduced in |AKMU IT] . Here we are 
only concerned with evolving preference lists (or permutations). In our model, we consider 
discrete time-steps. In each time-step, the algorithm is competing against nature as follows: 
the algorithm can query the input locally, nature lets the input evolve according to one or more 
evolution events. 

A query to the stable matching input is a triplet {z, re, re) G (Ax B‘^)U{B x A^) and the answer 
is whether 7r2(re) < 7r2(re). One evolution event consists of the following: pick 2 ; G A U R and 
i G [re — I] uniformly at random and swap re = Trf^{i) and re = Trf^{i 1) {i.e., set tTz{u) = i-\-l 
and t:z{v) = i). 
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While the algorithm can perform one query per time-step, a evolution events also occur, 
where a > 1 is some integer called the evolution rate. We further assume that a = 0(1), meaning 
that evolution events occur basically as often than the algorithm probes. We emphasize that 
the rate-limiting factor is the queries made by the algorithm. In particular, the algorithm may 
perform arbitrary (polynomial-time) computations in between time-steps. We are now ready to 
define our problem: 

Evolving Stable Matching (ESM): Given query access to an instance of the 
stable matching problem with evolution rate ol = 0(1), maintain a matching that 
minimizes the number of blocking pairs. 

2.3 Sorting evolving permntations 

The problem of sorting a single evolving permutation has been already addressed in [AKMUfT] . 
In this context, the evolution rate is still constant, but denotes the evolution rate of this single 
permutation. We will use the algorithm Quicksort of [AKMUfT] . It is simply the randomized 
version of quicksort which is shown to be robust with respect to an evolving input. The first 
lemma shows that the running time of quicksort is not affected by evolution events H 

Lemma 2.1 (Proposition 3 in |AKMIJfT| l. The running time of Quicksort is 0{nlogn) in 
expectation and with high probability, for any rate of evolution when the pairs to be swapped are 
chosen randomly. 

Second, Lemma 6 in |AKMUfT| states that Quicksort when run on an evolving permuta¬ 
tion TT, computes a permutation tt in which every element is approximately sorted. At time-step 
t, let vr* the denote the current permutation, and it* its approximation computed by the algo¬ 
rithm. 

Lemma 2.2 (Lemma 6 in [AKMU fl] 1 . Lett be the time-step of completion of Quicksort, then 
given an element u, the number of pairs (n, v) that the permutations vr* and tt* rank differently 
is O(logn) in expectation and with high probability. 

In our setting there are 2n evolving permutations over some set of n elements. Algorithm [T] 
simply sorts m (out of 2n) permutations, denoted by tti, ... ,TTm using Quicksort one after 
another. (We always invoke Algorithm [1] with either n or 2n permutations.) 


Algorithm 1 : Sequential sorting 


1: procedure SEQUENTlALSORT({7rj : j = 1,.. 
2; for j = 1 to m do 

3: TTj ^ QuiCKSORT(7rj) 

4: return {tTj : j = 1,..., m} 

., m}) > Only have query access to input 


Using Lemma 12.21 (Lemma 6 of |AKMUfT| l we can argue that Algorithm [T] maintains all 
permutations approximately sorted. While the evolving rate is still a = 0(1), there are now 2re 
evolving permutations, and the total number of evolution events is ol per time-step. 

Lemma 2.3. Let t be the time-step when Algorithm [I] terminates. Then, for m < 2n, given 
any element u and j G [m], the number of pairs {u,v) that the permutations iTj and vrj rank 
differently is O(logn) in expectation and with high probability. 

^We remark that Anagnostopoulos et al. | AKMUfT] use ‘whp’ to denote events that hold with probability 
1 — o(l), rather than the stronger notion we use in this paper. However, their proofs for the results used in our 
paper actually prove the stronger bounds. They have other results that do not satisfy the stronger notion and 
these are not used in our work. 
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Proof. Fix some j G [n]. Suppose that yrj was computed at time-step t' < t (the time-step when 
Quicksort for this particular list terminates). By Lemma 12.21 the statement holds for u at 
time-step t'. Due to Lemma l2.1l we have t — t' = 0{rPlogn) with high probability. During this 
time, the number of evolution steps that have swapped u with an adjacent element is O(logn) 
with high probability. This follows from a balls-and-bins experiment where we throw 0{n^ log n) 
balls (corresponding to the evolution steps) into m{n — 1) bins (corresponding to the adjacent 
pairs). It is known (see Exercise 3.1 in |MR95] 1 that in this particular case the number of balls 
in every bin is of the order of its mean with high probability. Therefore, during this time, at 
most O(logn) more elements may be swapped with u. □ 


2.4 Chernoff Bound with dependent variables 

We will require the following extension of the Chernoff bound. It follows from (the more general) 
Theorem 3.8 in |McD98| . 

Theorem 2.2. For i G [n], let Yi be a random variable over some set W and Xi be a Boolean 
random variable. For any y G 0?=! ^ ^ N? Ekiu) denote the event Yi = yi, ..., = 

yk. Suppose = l\Ek-i{y)] < p, for all y and k as above. Then, for any t > 0, 


P 


Xk>pn + t 


< exp 


3t^ \ 

6pn + 2t) 


3 Two simple cases 

In this section we present two simple arguments in two different settings. First, we consider 
how the original Gale-Shapley algorithm performs when run, without any modification, on lists 
produced by running quicksort on the evolving input. We present a simple analysis showing a 
bound of 0(n log n) on the number of blocking pairs, which is better than the trivial bound of 
order n^. Next, we analyze the Gale-Shapley algorithm when evolution events only occur on one 
side and the preference lists are uniformly random permutations; in particular the preference 
lists on the A side are chosen uniformly at random, and the preference lists on the B-side are 
subject to evolution events^ We present a simple analysis showing an O(logn) bound on the 
number of blocking pairs for this special case. 


3.1 A simple algorithm 

Our first algorithm ignores evolution of preference lists and runs the standard Gale-Shapley 
algorithm to produce a matching. More specifically, it first obtains the preferences lists for all 
2n agents using the Quicksort algorithm of |AKMU iT] [i.e., using Algorithm [T]) and then 
produces a matching by running the Gale-Shapley algorithm on these lists (ignoring the fact 
that these lists are only approximately correct). 

We show that this simple algorithm maintains a matching with at most 0(n log n) blocking 
pairs. Note that the number of blocking pairs is trivially at most n^. We further argue that 
improving the bound would require new ideas that either go around Lemma 12.31 (Lemma 6 of 
[AKMUfT] l or improve the analysis in a substantial way. 

Algorithm [2] runs in perpetuity. The matching M is maintained as the output until the new 
matching based on the newly sorted preference lists can be computed. 

® Note that after sufficiently many time-steps (though still polynomial) the evolution events ensure that all 
permutations are uniformly random. This follows from analyzing the mixing time of the corresponding Markov 
chain over permutations. See for example the book | LPWn9| . 
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Algorithm 2 : Simple dynamic stable matching 
1: while True do 

2: '■ Z ^ AU B} ^ SEQUENTiALSORT({7r^ : z G Au B}) > Calling Algorithm □ 

3: return Gale-Shapley matching M on the (approximately) sorted lists {tt^ : z G AU B} 


Theorem 3.1. For a sufficiently large constant cq and any time-step t > con^ logn, Algorithm\^ 
maintains a matching with 0(n log n) blocking pairs in expectation and with high probability. 


Proof. We consider the number of blocking pairs at time-step T < t when the current matching 
M was computed. In the following discussion, we use the following notation: for x G A,y G B, 
if M(x) = y, then M{y) = x (rather than M~^). At time-step T, for each z G Au B, define 
the indicator function Iz{w) to be 1 when M{z) w and w M{z) and 0 otherwise. (We 
don’t explicitly use superscripts on the preference lists vr as time-step T is fixed until specified 
otherwise.) By Lemma 12.31 in expectation and with high probability, 

= lire : rc M{z) and M{z) r(;}| = O(logn), (1) 

W 

If a pair {x,y) is blocking at time-step T, then x -<Try M{y) and y M{x). Assume x -<jry 
M{y) and y M[x). Since y FI(x}, x must have proposed to y at some point during the 
execution of the Gale-Shapley algorithm. By the properties of the Gale-Shapley algorithm, y 
should have been matched to an element of A with rank according to Hy at least as high as the 
rank of x in ity. This contradicts x -<jiy M{y). It follows that either M{y) x or M(x) -<jr^ y. 
Define U{x,y) to be 1 when (x,y) is blocking and 0 otherwise. We have argued that 

U{x,y) < Iffy) +Iy{x). 

By the union bound. Equation [1] holds for every z G Au B with high probability. Summing 
over all pairs {x, y) and applying the union bound again 

'^Uix^y) Iffy) + Iy{x) = Iffy)^ + Y (X] n), 

x^y x,y X y y x 

in expectation and with high probability. 

Next, let t' = t — T. We need to account for blocking pairs that may have arisen during 
these t' time-steps. First, we observe that by Lemma [2. II and union bound, with high probability 
t' = 0{n^ log n) (as otherwise another matching more up-to-date than the one at time-step T 
would be available). During the t' time-steps from T to t, evolution may create a blocking 
pair only if the swap decreases the rank of M{z) in vr^, for some z G Au B. Therefore, each 
step of the evolution introduces—independently—a new blocking pair with probability at most 
a/{n — 1). The expected number of blocking pairs introduced is therefore at most 0(nlogn) 
for a = 0(1) (and assuming t' = 0(n^ logn)). The result follows by a simple application of the 
Ghernoff and union bounds. □ 

Remark 3.1. We present a pair of instances to the stable matching problem with preference 
lists Tiz and tTz for z G A\J B, for which the conclusion of Lemma 12.31 is satisfied, i.e., for any z 
the number of pairs {iff) that are ordered differently in and tTz are O(logn). We then show 
that a matching that is stable with respect to tt as preference lists has D(nlogn) blocking pairs 
with respect to tt. Thus, it follows that using the Quicksort algorithm of Anagnostopoulos et 
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^-side lists vr B-side lists vr 


1 

3 2 1 4 5 6 7 


1 

6 7 1 5 4 3 2 

2 

4 3 2 5 6 7 1 


2 

7 1 2 6 5 4 3 

3 

5 4 3 6 7 1 2 


3 

1 2 3 7 6 5 4 

4 

6 5 4 7 1 2 3 


4 

2 3 4 1 7 6 5 

5 

7 6 5 1 2 3 4 


5 

3 4 5 2 1 7 6 

6 

1 7 6 2 3 4 5 


6 

4 5 6 3 2 1 7 

7 

2 1 7 3 4 5 6 


7 

5 6 7 4 3 2 1 


A-side lists vr i?-side lists it 


1 

1 2 3 4 5 6 7 


1 

1 7 6 5 4 3 2 

2 

2 3 4 5 6 7 1 


2 

2 1 7 6 5 4 3 

3 

3 4 5 6 7 1 2 


3 

3 2 1 7 6 5 4 

4 

4 5 6 7 1 2 3 


4 

4 3 2 1 7 6 5 

5 

5 6 7 1 2 3 4 


5 

5 4 3 2 1 7 6 

6 

6 7 1 2 3 4 5 


6 

6 5 4 3 2 1 7 

7 

7 1 2 3 4 5 6 


7 

7 6 5 4 3 2 1 


Figure 1: Instances vr and it demonstrating tightness of Algorithm [21 


al. [AKMUfTj and its analysis as a blackbox will not result in a stronger result than the one 
provided in Theorem 13.11 

First, we define the preference lists vt^ for z £ AU B. Let A = B = [n]. Then for x £ A, 
the preference list (ranking) is defined asx,x + l,...,n,l,2,...,x — 1. On the other hand 
for y £ B, the preference list (ranking) vr^ is defined as y,y — 1,... ,l,n,n — 1,... ,y + 1. The 
rankings vr^, z £ Au B, are now defined as follows: let k be some parameter, vr^ simply as the 
elements at rank 1 and k swapped. Figured] shows an example with n = 7 and k = 3. Clearly, 
when k = 0(logn), vr^ and vf^ satisfy the conclusion of Lemma 12.31 Yet, it is easy to see that 
M{i) = i is a stable matching for the preference lists vf^, z £ AU B, and for this matching with 
respect to the preference lists z £ AU B, every pair (i, j) with 0<j — i<k is a blocking 
pair. 


3.2 One-sided evolution 

In this section we analyze how the Gale-Shapley algorithm performs when the initial preference 
lists are random, but there is no evolution on the lists of the elements in A. Furthermore, we 
are going to assume that the algorithm knows each permutation in '■ x £ A}. We call this 
setting one-sided evolution. 


Algorithm 3 Gale-Shapley Algorithm 

1 

M 

> Initialize empty matching 

2 

for X £ A do 


3 

newjnatch £- False 


4 

p •(— X 


5 

while newjnatch = False do 


6 

y £- first as yet unproposed as per VTp 


7 

if M{y) not yet set then 


8 

M{p) £- y 


9 

M{y) 


10 

newjnatch = True 


11 

else if p -<x-y M{y) then 

t> y prefers p to M{y) 

12 

p' ^ M{y) 


13 

M{y) ^p 


14 

M{p) £- y 


15 

p p' 


16 

return M 



In this setting, the standard Gale-Shapley algorithm is implemented (basic pseudocode is 
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shown in Algorithm [3]) . Note that the only time the preference lists on the i?-side are used is in 
Line 11. Thus, it is only for these steps that we need to query the input (since the preference lists 
on the A-side are known to the algorithm). Thus, the number of queries made by the algorithm 
is bounded by the number of proposals. It was already observed that the number of proposals 
made in a random instance of the stable matching problem is O(nlogn). The actual algorithm 
keeps implementing the Gale-Shapley algorithm from scratch after completion. The matching 
from the previous completed run is used as the current matching. We prove the following result 
for the one-sided evolution setting. 

Theorem 3.2. For a sufficiently large constant cq and any time-step t > conlogn, the Gale- 
Shapley algorithm (repeatedly run and using the matching of the last completed run as the 
output) under one-sided evolution maintains a matching with at most O(logn) blocking pairs in 
expectation and with high probability. 


Proof. To prove the bound on the number of blocking pairs, call an evolution event on y’s list 
critical if it involves the then match of y. Suppose that after the algorithm terminates, y ^ B is 
involved in k blocking pairs (xi,y),..., (xfc,y). We observe that each one of the xi,... ,Xk was 
involved in at least one critical evolution step. To see this note that if {x, y) is blocking, then x 
proposed to y during the execution of the algorithm and got rejected subsequently (because tTx 
didn’t change and y M{x)). But since at the end of the execution it forms a blocking pair, 
it must ranked higher than M(y). This is only possible if x was swapped with the then match 
of y in some evolution event during the execution of the algorithm. 

Given this observation, we estimate the number of blocking pairs by estimating the number 
of critical evolution steps. Note that an evolution step is critical with probability at most 
2a/(n — 1) (at most 2n out of n(n — 1) pairs involve the matching). Let T be a random 
variable equal to the number of proposals before the algorithm outputs a matching and label 
the corresponding time-steps as 1,2, ...,r. For each step k, let be a Boolean random 
variable that is equal to 1 if at the time-step labeled k some evolution event was critical. 

First note that from a coupon-collecting argument as in Theorem 12.11 it follows that T = 
0(n log n) in expectation and with high probability. This is because for that argument the 
distribution of {vr^ : y € B} is irrelevant and {tt^ : x G A} being random permutations suffices. 
Therefore, we may fix an appropriately large constant C so that T < Cniogn with high 
probability and let m = Cniogn. As noted above, at any given step and given any information 
from the previous steps, an evolution step is critical with probability at most 2a/(n — 1); thus, 
by Theorem 12.21 


> 2Glogn 

k=l 


Oin '^), 


for sufficiently large C. We have 


T 

> 2Clogn 

k=l 


T 

< P [T > Gn log n] -|- P ^k > 2G log A (T < Cn log n) 

k=l 

m 

< P [T > Gn log n] -I- P > 2G log n . 

k=l 


By the observation at the beginning of this paragraph the claimed bound holds with high 
probability. The bound on the expectation follows by noting that there can be at most 
blocking pairs. 
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Finally, note that there will be at most 0(nlog n) time-steps before a new matching is 
computed. As in the final paragraph in the proof of Theorem 13.11 one can show that evolution 
cannot produce more than O(logn) blocking pairs in these many steps. □ 

4 General Case: Improved algorithm 

We now consider the general setting where the preference lists on both sides may be evolving. 
We present a modified version of the Gale-Shapley algorithm that takes advantage of Lemma [2. 3 1 
(Lemma 6 of [AKMUlT] i and maintains a stable matching with at most 0((logn)^) blocking 
pairs. The algorithm consists of two separate processes that run in an interleaved fashion: 
the sorting process on even time-steps and the matching process on odd ones. The sorting 
process is basically a call to SEQUENTiALSORxdvTa; | x G A}) that produces approximately 
sorted preference lists on the A side, {itj, | x G A}. The algorithm runs in perpetuity, in the 
sense that as soon as it terminates it restarts, though the copies {^x \ x G A} from the previous 
execution are retained to be used by the stable matching process. Initially, the tTx are set to be 
random permutations, thus, for the first O(re^logn) steps, until one run of the sorting process 
is complete, the matching output by the algorithm will be garbage. 

Algorithm 4 : Interleaving Sorting and Matching 
1: for t = 1, 2,... do 
2: if t is Even then 

3: Perform query for Algorithm [1] 

4: else if t is Odd then 

5: Perform query for Algorithm [5] 


We note that what is counted here is only time-steps; making queries is the bottleneck. 
Additional computations required by the algorithm can be performed in between time-steps at 
no additional cost. 

The sorting process performs queries only during even steps and its purpose it to keep the 
preference list of each x G A approximately sorted, where by approximately sorted we meant 
that the conclusion of Lemma 12.31 holds. 

The matching process performs queries during odd steps. Our stable matching algorithm, 
which is a variant of the Gale-Shapley algorithm, is presented in Algorithm [5l Note that the 
{tTx I X G A} used in Algorithm [5] are static and are the output of the latest completed run of 
Algorithm [TJ However, the comparisons made are all on dynamic data. The difference from the 
standard Gale-Shapley algorithm is that whenever some x G A is about to make a proposal, 
first the best y ^ B among the O(logn) highest ranked as per the ranking tTx that have not 
yet rejected x is found. Note however, that the best is with respect to the dynamic (current) 
preference list tTx (otherwise, it would be trivial since tTx is static). This operation is basically 
the algorithm to hnd the minimum element, which can be implemented in O(logn) time using 
only comparison queries (see Section 3 in |AKMU iT] 1. We don’t need to use any particular 
result regarding finding the minimum element; instead, we incorporate the errors that may 
have occurred while finding the minimum due to the dynamic nature of the data, into our 
stable matching analysis directly. 

We first describe the high-level idea of the proof. The sorting process needs 0(n^ log n) 
comparisons with high probability. The approximations {^x '■ x G A} of {tTx '■ x G A} that are 
being computed are used by the modified Gale-Shapley algorithm. By Lemma 12.31 we are able 
to claim that for any element u in the preference list of any x G A, the number of pairs (tt, v) 
that are ordered differently in tTx and tTx are O(logn). Therefore, when x is about to propose it 
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Algorithm 5 : Modified Gale-Shapley Algorithm 
1: M^0 
2: for X G A do 
3: newjnatch •(— FALSE 

4: P ^ X 

5: while newjnatch = False do 

6: 5 ■(—Clogn highest-ranked, not-yet-proposed elements in S per Tfj; 

7: y best(5) t> Best with respect to dynamic 

8: if M(y) not yet set then 

9: M{p) y 

10: M{y) ^ p 

11: newjnatch •<— True 

12: else if p M{y) then 

13: p' ^ M{y) 

14: M{y) ^ p 

15: M{p) y 

16: p ^ p' 

17: return M 


suffices to look among 0(log n) elements in ^x to find the y to which the proposal will be made. 
Since the matching process is expected to make 0(n log n) proposals, it is expected to require 
0(n(logn)^) comparisons. It turns out that, during these steps, evolution creates a blocking 
pair with probability at most a/n. Therefore we expect 0((logn)^) blocking pairs. 

We now provide the details of the proof. In order to bound the number of blocking pairs, 
it is crucial that during the matching process not too many queries are made, or alternatively 
that not too many proposals are made. We therefore need an analog of Theorem 12.11 To apply 
the coupon-collecting argument from the proof of that theorem we prove the following lemma. 

Lemma 4.1. Provided tXx was chosen uniformly at random at time-step 0 and only comparison 
queries are made, the element y chosen at line 7 of Algorithmic is a random element from the 
subset of B to which x has not by that point made any proposals. 

Proof. The proof of the lemma relies on the fact that the dynamic quicksort algorithm used to 
obtain ^x for x £ A and the procedure used to find the best element in line 7 of Algorithm [5] 
only use comparison queries. 

Let TTx be the preference list of x before the first comparison is queried. Fix an arbitrary 
sequence of evolution steps that occurs during the computation of y. Suppose that given these 
choices of nature, y = 7rx{k). Then, given the same evolution steps, for any other permutation 
y = Tr'xik). Since tTx is a random permutation, y = TTx{k) is a random element of era,. □ 

Remark 4.1. We remark that the requirement on the implementation of Quicksort and 
best using only comparison queries is necessary and the lemma does not hold for an arbitrary 
algorithm. 

Lemma 4.2. The number of proposals during one execution of the matching process is 0{n log n) 
in expectation and with high probability. 

Proof. Suppose x proposes to y and at that point k elements of B are unmatched. Note that 
it must be the case that x has not proposed to any of these elements (otherwise, they would 
not be currently unmatched). Thus, by the previous lemma, each of these k elements receives 
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a proposal with probability at least 1/n. The stated bound follows from the analysis of coupon 
collector’s problem as in Theorem 12.11 □ 


As in the one-sided setting, the analysis will rely on estimating the occurrence of a specific 
kind of critical evolution steps. In the present case the definition of a critical evolution step is 
a little more involved than its one-sided counterpart. 

Definition 4.1. An evolution event on the preference list vr^ of z a Au B is critical if one of 
the following holds: 

1. An evolution event involves a swap of the then mateh of z, M{z). 

2. If z £ A, the evolution event involves swapping the then best element as per vr^ to which 
z has not yet proposed. 

The following claim establishes the link between the critical evolution steps and the number 
of blocking pairs. 

Claim 4.1. Assume that for the duration of one run of the matching process, the preference 
lists {t^z I G Au B} satisfy the conditions of Lemma \2.A. and suppose that {x,y) is a blocking 
pair with respect to the returned matching. Then there was a critical evolution event on the 
preference list of at least one of x and y during the execution of the matching process. 

Proof. First consider the case that x proposed to y during the execution of the matching process. 
It follows that y rejected x at some point in favor of some other element. When x was rejected, 
the then M{y) satisfies M{y) x. Subsequently, M{y) may change but has to become better, 
unless there was a swap that involves the then M{y), which is a critical event on tt^. Since, 
we know that in the end x -<Try x' , where x' is the final match of y, there must have been 
some evolution event where x was swapped with the then match of y. Thus, by part 1 of 
Definition SH a critical evolution event occurred on iTy. 

On the other hand, suppose x never proposed to y, and let y' be the final match of x. 
Suppose that when x proposed to y', y y'\ it follows that best on line 7 of Algorithm [5] 
failed to return the best element to which x had not yet proposed. Since, we are assuming that 
TTx is a sufficient approximation of TTa,, it must be because the actual best element was swapped at 
least once while best was being executed. Thus, by part 2 of Definition 14.11 a critical evolution 
event occurred on tTx. Finally, if when x proposed to y' , it was the case that y' ^7r„, y, but (x, y) 
is a blocking pair, it must be that y' was involved in a swap subsequently leading to a critical 
event involving x’s the then match. □ 

Theorem 4.1. Provided the initial preference lists are drawn randomly]^ for all z & AVJ B, 
for a sufficiently large constant cq and any time-step t > log n, Algorithm 0 maintains a 
matching with at most 0((logn)^) blocking pairs in expectation and with high probability. 


Proof. As a result of Claim 14.11 we can estimate the number of blocking pairs by estimating 
the number of critical evolution steps. Let T be a random variable equal to the number of 
queries of Algorithm [5] before it outputs a matching and label the corresponding time-steps as 
1,2,... ,T. For each step k, let be a Boolean random variable that is equal to 1 if at the 
time-step labeled k the evolution was critical. Furthermore, denote by £ the event that during 
these T time-steps the lists were approximately sorted. By Lemma 12.31 the event £ occurs with 
high probability. 

^This is not actually required, since after sufficiently long (though still polynomial) time, all the preference 
lists will be close to random due to a mixing time argument on the set of permutations. 
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By Lemma 14.21 it follows that T = 0(n(logn)^) in expectation and with high probability, 
since we are wasting O(logn) queries per proposal. Therefore, we may fix an appropriately 
large constant C so that T < C'n(logn)^ with high probability and let m = C'n(logn)^. Note 
that—given any history of evolution steps—an evolution step is critical with probability at most 
0{a/n), since for each z a Au B there is a constant number of elements that evolution has to 
swap in order to be critical. Thus, by Theorem 12.21 


p[^Xfc>2C(logn)2 

k=l 


0{n ^), 


for sufficiently large C. We have 


1 

Xk > 2C(log nf < P [T > Cn(log nf] + F[£] 


k=l 


T 

+ > 2C'(logn)^^ a(T < C'n(logn)^) A £ 

k=l 

m 

< P[T > Cn(logn)2] +P[,f] +pJ^Xfc > 2C{lognf . 

k=l 

It follows that the claimed bound holds with high probability. The bound on the expectation 
follows by noting that there can be at most blocking pairs. □ 
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